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Abstract 

The scaling properties encompass in a simple analysis many of the volatility charac- 
teristics of financial markets. That is why we use them to probe the different degree 
of markets development. We empirically study the scaling properties of daily For- 
eign Exchange rates, Stock Market indices and fixed income instruments by using 
the generalized Hurst approach. We show that the scaling exponents are associated 
with characteristics of the specific markets and can be used to differentiate mar- 
kets in their stage of development. The robustness of the results is tested by both 
Monte-Carlo studies and a computation of the scaling in the frequency-domain. 
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1 Introduction 



In a recent book [Dacorogna et al., 2001a], the hypothesis of heterogeneous 
market agents was developed and backed by empirical evidences. According 
to this view, the agents are essentially distinguished by the frequency at which 
they operate in the market. The scahng analysis, which looks at the volatihty 
of returns measured at different time intervals, is a parsimonious way of as- 
sessing the relative impact of these heterogeneous agents on price movements. 
Viewing the market efficiency as the result of the interaction of these agents 
[Dacorogna et al., 2001b], brings naturally to think that it is the presence of 
many different agents that would characterize a mature market, while the ab- 
sence of some type of agents should be a feature of less developed markets. 
Such a fact should then reflect in the measured scaling exponents. The study 
of the scaling behaviors must therefore be an ideal candidate to characterize 
markets. 

For institutional investors, a correct assessment of markets is very important to 
determine the optimal investment strategy. It is common practice to replicate 
an index when investing in well developed and liquid markets. Such a strategy 
minimizes the costs and allow the investor to fully profit from the positive 
developments of the economy while controlling the risk through the long ex- 
perience and the high liquidity of these markets. When it comes to emerging 
markets, it is also clear that the stock indices do not fully represent the un- 
derlying economies. Despite its higher costs, an active management strategy 
is required to control the risks and fully benefit from the oportunities offered 
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by these markets. The differentiation between markets is clear for the extreme 
cases: New York stock exchange and the Brazihan or Russian stock exchange. 
The problem lies for all those in between: Hungary, Mexico, Singapore and 
others. For those markets a way to clarify the issue will help decide on the 

best way to invest assets. 

The purpose of this article is to report on the identification of a strong rela- 
tion between the scaling exponent and the development stage of the market. 
This conclusion is backed by a wide and unique empirical analysis of several 
financial markets (32 Stock market indices, 29 Foreign exchange rates and 28 
fixed income instruments) at different development stage: mature and liquid 
markets, emerging and less liquid markets. Furthermore, the robustness and 
the reliability of the method is extensively tested through several numerical 
tests, Monte Carlo simulations with a variety of random generators and the 
comparison with results obtained from a frequency domain computation of 
related exponents. 

The scaling concept has its origin in physics but is increasingly applied out- 
side its traditional domain [Miiller et al., 1990] [Dacorogna et al., 2001a]. In 
the recent years, its application to financial markets, initiated by Mandelbrot 
in the 1960 [Mandelbrot, 1963, Mandelbrot, 1997], has largely increased also 
in consequence of the abundance of available data [Miiller et al., 1990]. Two 
types of scaling behaviors are studied in the finance literature: 

(1) The behavior of some forms of volatility measure (variance of returns, 
absolute value of returns) as a function of the time interval on which the 
returns are measured. (This study will lead to the estimation of a scaling 
exponent related to the Hurst exponent.) 
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(2) The behavior of the tails of the distribution of returns as a function of the 
size of the movement but keeping the time interval of the returns con- 
stant. (This will lead to the estimation of the tail index of the distribution 
[Dacorogna et al., 2001a].) 

Although related, these two analysis lead to different quantities and should 
not be confused as it is often the case in the literature as can be seen in the 
papers and debate published in the November 2001 issue of Quantitative Fi- 
nance [LeBaron, 2001] [Lux, 2001] [Mandelbrot, 2001]. For more explanations 
about this and the relation between the two quantities, the reader is referred 
to the excellent paper by [Groenendijk et al., 1998]. In this study, we are in- 
terested in the first type of analysis. Until now, most of the work has concen- 
trated in studies of particular markets: Foreign Exchange [Miiller et al., 1990] 
[Dacorogna et al, 2001a] [Corsi et al., 2001], US Stock Market (Dow Jones) 
[Mantegna and Stanley, 1995] or Fixed Income [Ballocchi et al., 1999]. These 
studies showed that empirical scaling laws hold in all these markets and for a 
large range of frequencies: from few minutes to few months. 

Recently, a controversy has erupted between [LeBaron, 2001] on one side 
and [Mandelbrot, 2001, Stanley and Plerou, 2001] on the other side with some- 
where in the middle [Lux, 2001] to know if the processes that describe financial 
data are truly scaling or simply an artifact of the data. Moreover, these papers 
propose new scaling models or empirical analysis that better describe empiri- 
cal evidences and one could add to these [Bouchaud et al., 2000]. It should be 
however noted that - as underhned by Stanley et al. [Stanley et al., 1996] - in 
statistical physics, when a large number of microscopic elements interact with- 
out characteristic scale, universal macroscopic scaling laws may be obtained 
independently of the microscopic details. 
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Here we address the question of the scahng properties of financial time series 
from another angle. We are not interested in fitting a new model but want to 
gather empirical evidences by analyzing daily data (described in Section 2). 
With the same methodology, we study very developed as well as emerging 
markets in order to sec if the scaling properties differ between the two and if 
they can serve to characterize and measure the development of the market. 
Here the scaling law is not used to conclude anything on the theoretical process 
but to the contrary we use it as a "stylized fact" that any theoretical model 
should also reproduce. Our purpose is to show how a relatively simple statistics 
gives us indications on the market characteristics, very much along the lines of 
the review paper by Brock [Brock, 1999] . In Section 3, we recall the theoretical 
framework and in Section 4 we introduce the generalized Hurst exponents 
analysis. The metodology is described in Section 5. In Section 6 the generalized 
Hurst exponents results and their temporal stability check are presented. In 
Section 7 we compute the scaling exponents in the frequency domain and we 
compare the scaling spectral exponents and the Hurst exponents. In Section 8 
a Monte Carlo simulation is presented. Finally some conclusions are given in 
Section 9. 

2 Data Description and Studied Markets 

We study several financial markets which are at different development stage: 
mature and liquid markets, emerging and less liquid markets. Moreover, we 
choose markets that deal with different instruments: equities, foreign exchange 
rates, fixed income futures. In particular, the data that we analyze are: twenty 
nine Foreign Exchange rates (FX) (see Table 1), thirty two Stock Market 
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indices (SM) (see Table 2), Treasury rates corresponding to twelve different 
maturity dates (TR) (see Table 3) and Eurodollar rates having maturity dates 
ranging from 3 months to 4 years (ER) (see Table 4). Hereafter we give a brief 
description of the time-series studied in this paper. 

FX: The Foreign Exchange rates (Table 1) are 29 daily spot rates of major 
currencies against the U.S. dollar. The time series that we study go from 
1990 to 2001 and 1993 to 2001. These rates have been certified by the Federal 
Reserve Bank of New York for customs purposes. The data are noon buying 
rates in New York for cable transfers payable in the hsted currencies. These 
rates are also those required by the Securities and Exchange Commission 
(SEC) for the integrated disclosure system for foreign private issuers. The 
information is based on data collected by the Federal Reserve Bank of New 
York from a sample of market participants. 

SM: The Stock Market indices (reported in Table 2) are 32 of the major 
indices of both very developed markets like the US or European markets 
and emerging markets. These daily time series range from 1990 or 1993 to 
2001. 

TR: The Treasury rates (Table 3) are daily time series going from 1990 to 
2001. The yields on Treasury securities at 'constant maturity' are interpo- 
lated by the U.S. Treasury from the daily yield curve. This curve, which 
relates the yield on a security to its time to maturity, is based on the clos- 
ing market bid yields on actively traded Treasury securities in the over- 
the-counter market. These market yields are calculated from composites of 
quotations obtained by the FD Bank of New York. The constant maturity 
yield values are read from the yield curve at fixed maturities, currently 3 
and 6 months and 1, 2, 3, 5, 7, 10, and 30 years. The Treasury bill rates 
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are based on quotes at the official close of the U. S. Government securities 
market for each business day. They have maturities of 3 and 6 months and 
1 year. 

ER: The Eurodollar interbank interest rates (Table 4) are bid rates with 
different maturity dates and they are daily data in the time period 1990- 
1996 [Di Matteo and Aste, 2002]. 

3 Theoretical Framework and background 

The scaling properties in time series have been studied in the literature by 
means of several techniques. For the interested reader we mention here some 
of them such as the seminal work [Hurst, 1951] on rescaled range statis- 
tical analysis R/ S with its complement [Hurst et al., 1965] and the modi- 
fied R/S analysis of [Lo, 1991], the multiaffine analysis [Peng et al., 1994], 
the dctrcndcd fluctuation analysis (DFA) [Ausloos, 2000], the pcriodogram 
regression (GPH method) [Geweke and Porter- Hudak, 1983] , the (m, /c)-Zipf 
method [Zipf, 1949], the moving- average analysis technique [EUinger, 1971], 
the Average Wavelet Coefficient Method in [Percival and Walden, 2000] and 
in [Gengay et al., 2001], the ARFIMA estimation by exact maximum likeli- 
hood (ML) [Sowell, 1992] and connection to multi-fractal/multi-affine anal- 
ysis (the q order height-height correlation) have been made in various pa- 
pers like [Ivanova and Ausloos, 1999] . In the financial and economic litera- 
ture, many are the proposed and used estimators for the investigation of the 
scaling properties. To our knowledge it doesn't exit one whose performance 
has no defeciencies. The use of each of the above mentioned estimators can 
be subject to both advantages and disadvantages. For instance, simple tradi- 
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tional estimators can be seriously biased. On the other hand, asymptotically 
unbiased estimators derived from Gaussian ML estimation are available, but 
these are parametric methods which require a parameterized family of model 
processes to be chosen a priori, and which cannot be implemented exactly in 
practice for large data sets due to high computational complexity and memory 
requirements [Phillips, 1999a,Phillips, 1999b, Phillips, 2001] . Analytic approx- 
imations have been suggested (Whittle estimator) but in most of the cases 
(see [Beran, 1994]), computational difficulties remain, motivating a further 
approximation: the discretization of the frequency-domain integration. Even 
with all these approximations the Whittle estimator remains with a signifi- 
cantly high overall computational cost and problems of convergence to local 
minima rather than to the absolute minimum may be also encountered. 

The rescaled range statistical analysis [R/ S analysis) was first introduced 
by Hurst to describe the long-term dependence of water levels in rivers and 
reservoirs. It provides a sensitive method for revealing long-run correlations in 
random processes. This analysis can distinguish random time series from cor- 
related time series and gives a measure of a signal "roughness" . What mainly 
makes the Hurst analysis appealing is that all these information about a com- 
plex signal are contained in one parameter only: the Hurst exponent. How- 
ever, the original Hurst R/S approach has problems in the presence of short 
memory, heteroskedasticity, multiple scale behaviors. This has been largely 
discussed in the literature (see for instance [Lo, 1991,Teverovsky et al., 1999]) 
and several alternative approaches have been proposed. The fact that the 
range relies on maxima and minima makes also the method error-prone to 
any outlier. Lo [Lo, 1991] suggested a modified version of the R/S analy- 
sis that can detect long-term memory in the presence of short-term depen- 
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dence [Moody and Wu, 1996]. The modified R/S statistic differs from the 
classical R/S statistic only in its denominator, adding some weights and co- 
variance estimators to the standard deviation suggested by Newey and West 
[Newey and West, 1987], and a truncation lag, q. In the modified R/S, a prob- 
lem is choosing the truncation lag q. Andrews [Andrews, 1991] showed that 
when q becomes large relative to the sample size N, the finite-sample distri- 
bution of the estimator can be radically different from its asymptotic limit. 
However, the value chosen for q must not be too small, since the autocorrela- 
tion beyond lag q may be substantial and should be included in the weighted 
sum. The truncation lag thus must be chosen with some consideration of the 
data at hand. 

In this paper we use a different and alternative method: the generalized Hurst 
exponent method. We choose this type of analysis precisely because it com- 
bines the sensitivity to any type of dependence in the data to a computation- 
ally straight forward and simple algorithm. The main aim of this paper is to 
give an estimation tool from an empirical analysis which provides a natural, 
unbiased, statistically and computationally efficient, estimator of the gener- 
alized Hurst exponents. This method, described in the following section, is 
first of all a tool which studies the scahng properties of the data directly 
via the computation of the g-order moments of the distribution of the incre- 
ments. The g-order moments are much less sensitive to the outliers than the 
maxima/minima and different exponents q are associated with different char- 
acterizations of the multiscaling complexity of the signal. In the following we 
show that this method is robust and it captures very well the scaling features 
of financial fiuctuations. We show that through the use of a relatively simple 
statistics we give a wide view of the scaling behaviour across different markets. 
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4 Generalized Hurst exponent 



The Hurst analysis examines if some statistical properties of time series X{t) 
(with t—u, 2v, kv, T) scale with the observation-period (T) and the 
time-resolution {y) . Such a scaling is characterized by an exponent H which is 
commonly associated with the long-term statistical dependence of the signal. 
A generalization of the approach proposed by Hurst should therefore be associ- 
ated with the scaling behavior of statistically significant variables constructed 
from the time series. To this purpose we analyze the g-order moments of the 
distribution of the increments [Mandelbrot, 1997,Barabasi and Vicsek, 1991] 
which is a good characterization of the statistical evolution of a stochastic 
variable X{t): 

where the time- interval r can vary between v and Tmax- (Note that, for q — 2, 
the Kq{T) is proportional to the autocorrelation function: a(T) = {X{t + r)X(t)).) 

The generalized Hurst exponent H(q) ^ can be defined from the scaling be- 
havior of Kq{T) [Barabasi and Vicsek, 1991], which can be assumed to follow 
the relation 



K,{r) ~ (-) . (2) 



This assumption flows naturally from the result of [Groenendijk et al., 1998] 
and has been carefully checked to be hold for the financial time series studied in 



^ We use H without parenthesis as the original Hurst exponent and H{q) as the 
generalized Hurst exponent. 
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this paper. For instance, in Fig. 1, the scahng behaviour of Kq{T) in agreement 
with Eq. 2 is shown in the time period from 1990 to 2001 for Nikkei 225. Each 
curve corresponds to different fixed values of q ranging from q = 1 to q — 3, 
whereas r varies from 1 day to 19 days. 

Within this framework, we can distinguish between two kind of processes: (i) 
a process where H[q) = H, constant independent of q; (ii) a process with 
H{q) not constant. The first case is characteristic of uni-scaling or uni-fractal 
processes and its scaling behavior is determined from a unique constant H 
that coincides with the Hurst exponent. This is for instance the case for self- 
affine processes where qH{q) is linear {H{q) — H) and fully determined by 
its index H. (Recall that, a transformation is called affine when it scales time 
and distance by different factors, while a behavior that reproduces itself under 
affine transformation is called self-affine [Mandelbrot, 1997]. A time-dependent 
self-affine function X{t) has fluctuations on different time scales that can be 
rescaled so that the original signal is statistically equivalent to its rescaled 
version \~^X{\t) for any positive A, i.e. X{t) ~ X~^X{\t). Brownian mo- 
tion is self-affine by nature.) In the second case, when H[q) depends on q, 
the process is commonly called multi-scaling (or multi-fractal) and different 
exponents characterize the scaling of different g-moments of the distribution. 

For some values of q, the exponents are associated with special features. For 
instance, when g = 1, H{1) describes the scaling behavior of the absolute 
values of the increments. The value of this exponent is expected to be closely 
related to the original Hurst exponent, H, that is indeed associated with the 
scaling of the absolute spread in the increments. The exponent at g 2, is 
associated with the scaling of the autocorrelation function and is related to 
the power spectrum [Flandrin, 1989]. A special case is associated with the 
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value of q — q* at which q*H{q*) = 1. At this value of q, the moment Kg*{T) 
scales linearly in r [Mandelbrot, 1997]. Since qH{q) is in general a monotonic 
growing function of g, we have that all the moments Hg^r) with q < q* will 
scale slower than r, whereas all the moments with q > q* will scale faster 
than r. The point q* is therefore a threshold value. In this paper wc focal- 
ize the attention on the case q — 1 and 2. Clearly in the uni- fractal case 
H{1) — H{2) — H{q*). Their values will be equal to 1/2 for the Brownian 
motion and they would be equal to H ^ 0.5 for the fractional Brownian mo- 
tion. However, for more complex processes, these coefficients do not in general 
coincide. We thus see that the non-linearity of the empirical function qH{q) is 
a solid argument against Brownian, fractional Brownian, Levy, and fractional 
Levy models, which are all additive models, therefore giving for qH{q) straight 
lines or portions of straight lines. The curves for qH{q) vs. q are reported in 
Fig. 2 for some of the data. One can observe that, for all these time series, 
qH{q) is not linear in q but slightly bending below the linear trend. The same 
behavior holds for the other data. This is a sign of deviation from Brownian, 
fractional Brownian, Levy, and fractional Levy models, as already seen in FX 
rates [Miiller et al., 1990]. 

5 Methodology and Preparation of the Data 

Let us here recall that the theoretical framework we presented in the previous 
section is based on the assumption that the process has the scaling property 
described in Eq. 2. Moreover, we have implicitly assumed that the scaling 
properties associated with a given time series stay unchanged across the ob- 
servation time window T. On the other hand, it is well known that financial 
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time series show evidences of variation of their statistical properties with time, 
and show dependencies on the observation time window T. The simplest case 
which shows such a dependence is the presence of a linear drift {r]t) added 
to a stochastic variable {X{t) — X{t) + r}t) with X{t) satisfying Eq. 2 and 
the above mentioned properties of stability within the time window. Clearly, 
the scaling analysis described in the previous section must be applied to the 
stochastic component X{t) of the process. This means that we must subtract 
the drift rjt from the variable X{t). To this end one can evaluate 7] from the 
following relation: 

(X(i + T)-X(i))=77T . (3) 



Other more complex deviations from the stationary behavior might be present 
in the financial data that we analyze. In this context, the subtraction of the 
linear drift can be viewed as a first approximation. 

Our empirical analysis is performed on the daily time series TR, ER, FX and 
SM (described in Section 2) which span typically over periods between 1000 
and 3000 days. In particular, we analyze the time series themselves for the TR 
and ER, whereas we compute the returns from the logarithmic price X{t) — 
ln(P(t)) for FX and SM. Moreover, all of these variables arc 'dctrended' by 
eliminating the linear drift (if there is one) as described in (Eq. 3). 

We compute the g-order moments Kq{T) (defined in Eq. 1) of the 'detrended' 
variables and their logarithms with r in the range between v — 1 day and 
Tjnax days. In order to test the robustness of our empirical approach, for each 
series we analyze the scaling properties varying T^ax between 5 and 19 days. 
We compute the 99% confidence intervals of all the exponents using differ- 
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ent Tmax values ^ . The resulting exponents computed using different Tmax are 
stable in their values within a range of 10%. We then verify that the scal- 
ing behavior given in Eq. 2 is well followed (see Fig. 1) and we compute the 
associated generalized Hurst exponent H{q) whose values are given in the fol- 
lowing section. We also tested the influence of the detrending (through Eq. 3) 
calculating the generalized Hurst exponent both for the detrended and the 
non-detrending time series. The results are in all cases comparable within the 
standard deviations calculated varying T^ax- 

6 Results 

6. 1 Computation of the generalized Hurst exponent 

In this section we report and discuss the results for the scaling exponents H{q) 
computed for g = 1 and q — 2. These exponents H(l) and H(2) for all the as- 
sets and different markets (presented in Section 2) are reported in Figs. 3 and 4 
respectively. Figures 3 (a) and 4 (a) refer to the Treasury and Eurodollar rates 
in the time period from 1990 to 1996. Whereas Figures 3 (b) and 4 (b) are 
relative to the Stock Market indices and Foreign Exchange rates in the time 
period reported in Tables 1 and 2. The data points are the average values 
of H{1) and H{2) computed from a set of values corresponding to different 
Tmax (between 5 and 19 days) and the error bars are their standard deviations. 
The generalized Hurst exponents are computed through a linear least squares 
fitting. We have computed the standard deviations for the two linear fit coef- 

^ By using a Matlab routine, namely, normfit that computes parameter estimates 
and confidence intervals for normal data. 
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ficients and the correlation coefficient. It results that the standard deviations 
from the linear fitting are below or equal to the reported standard deviations 
values computed varying Tmax- The correlation coefficient is never lower that 
0.99. 

Let us first notice that, for fixed income instruments (Figs. 3(a) and 4(a)), 
H{2) is close to 0.5 while H[l) is rather systematically above 0.5 (with the 
3 months Eurodollar rate that shows a more pronounced deviation because 
it is directly infiuenced by the actions of central banks). On the other hand, 
as far as Stock markets are concerned, we find that the generalized Hurst 
exponents H{1), H{2) show remarkable differences between developed and 
emerging markets. In particular, the values of H{1), plotted in Fig. 3(b), 
present a differentiation across 0.5 with high values of H{1) associated with 
the emerging markets and low values of H{1) associated with developed ones. 
In Fig. 3(b) the ordering of the stock markets from left to right is chosen 
in ascending order of H{1). One can see that such a ordering corresponds 
very much to the order one would intuitively give in terms of maturity of the 
markets. Moreover, we can see from Fig. 4(b) that the different assets can be 
classified into three different categories: 

(1) First those that have an exponent H{2) > 0.5 which includes all indices 
of the emerging markets and the BCI 30 (Italy), IBEX 35 (Spain) and 
the Hang Seng (Hong Kong). 

(2) A second category concerns the data exhibiting H{2) ~ 0.5 (within the 
error bars). This category includes: FTSE 100 (UK), AEX (Nether- 
lands), DAX (Germany), Swiss Market (Switzerland), Top 30 Capital 
(New Zealand), Tel Aviv 25 (Israel), Seoul Composite (South Korea) and 
Toronto SE 100 (Canada). 
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(3) A third category is associated with H{2) < 0.5 and includes the following 
data: Nasdaq 100 (US), SkPbOO (US), Nikkei 225 (Japan), Dow Jones 
Industrial Average (US), CAC 40 (France) and All Ordinaries (Australia). 

We find therefore that all the emerging markets have H{2) > 0.5 whereas all 
the well developed have H{2) < 0.5. This simple classification is not achieved 
by other means. One could, for instance use the Sharpe Ratio [Sharpc, 1994] 
we have tried it but it does not achieve such a clear cut categorization. This 
ratio requires a benchmark risk free return that is not always available for 
emerging markets. We have tried the classification obtained from a simple 
ratio of the average returns of their standard deviations but the ordering is 
not conclusive. 

We find that the Foreign Exchange rates show H{1) > 0.5 quite system- 
atically. This is consistent with previous results computed with high fre- 
quency data [Miiller et al., 1990], although the values here are shghtly lower. 
An exception with pronounced H{1) < 0.5 is the HKD/USD (Hong Kong) 
(Fig. 3 (b)). This FX rate is, or has been, at one point pegged to the USD, that 
is why its exponent differs from the others. Whereas in the class H(l) ~ 0.5 
we have: IT L /USD (Italy), PHP /USD (Philippines), AUD/USD (Aus- 
tralia), NZD/USD (New Zealand), ILS/USD (Israel), CAD /USD (Canada), 
SGD/USD (Singapore), NLG/USD (Netherlands) and JPY/USD (Japan). 
On the other hand, the values of H{2) (Fig. 4 (b)) show a much larger ten- 
dency to be < 0.5 with some stronger deviations such as: HKD/USD (Hong 
Kong), PHP /USD (Philippines), KRW/U SD{^outh Korea), PEN /USD 
(Peru) and TRL/USD (Turkey). Whereas values of H{2) > 0.5 are found in: 
GBP/USD (United Kingdom), PESO /USD (Mexico), INR/USD (India), 
IDR/USD (Indonesia), TWD/USD (Taiwan) and BRA/USD (Brazil). 
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6.2 Checking the Temporal and Numerical Stability of the Results 

In order to check the temporal stabihty of the results, these analysis are per- 
formed also over different time periods and the values of the exponents H{1) 
and H{2) are reported in Table 6 for the time period from 1997 to 2001 for For- 
eign Exchange rates, Stock Market indices and Treasury rates. These results 
should be compared with those obtained on the whole time period (shown in 
Table 2) and to time periods of 250 days. Moreover, we tested the numeri- 
cal robustness of our results by using the Jackknife method [Kunsch, 1989] 
which consists of taking out randomly 1/10 of the sample and iterates the 
procedure 10 times (every time taking out data which were not taken out 
in the previous runs). On one hand, we observe (see Fig. 5) that the gener- 
alized Hurst exponents computed on these Jackknife-reduced time series are 
very close to those computed on the entire series with deviations inside the 
errors estimated by varying Tmax (as described in section 5). This indicates 
a strong numerical stability. On the other hand, the analysis on sub-periods 
of 250 days shows fluctuations that are larger than the previous estimated 
errors (and larger than the variations with the Jackknife method) indicat- 
ing therefore that there are significant changes in the market behaviors over 
different time periods (Fig. 5 (a)). This phenomenon was also detected in 
[Dacorogna et al., 2001a] when studying Exchange rates that were part of the 
European Monetary System. It seems that H{1) is particularly sensitive to in- 
stitutional changes in the market. The scaling exponents cannot be assumed to 
be constant over time if a market is experiencing major institutional changes. 
Nevertheless, well developed markets have values of H{2) that are on aver- 
age smaller than the emerging ones and the weakest markets have oscillation 
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bands that stay above 0.5 whereas the strongest have oscillation bands that 
contain 0.5. 

Numerical stability gives us confidence in our method for determining the 
exponent and temporal variability is a sign that the exponents are sensitive 
to institutional changes in the market reinforcing our idea to use them as 
indicators of the maturity of the market. 

7 Scaling exponents in the frequency domain 

7.1 Spectral analysis 

In order to empirically investigate the statistical properties of the time series 
in the frequency domain we perform a spectral analysis computing the power 
spectral density (PSD) [Kay and Marple, 1981] by using the periodogram ap- 
proach, that is currently one of the most popular and computationally efficient 
PSD estimator. This is a sensitive way to estimate the limits of the scaling 
regime of the data increments. The results for some SM data in the time pe- 
riods 1997-2001 are shown in Fig. 6. For SM we compute the power spectra 
of the logarithm of these time series. As one can see the power spectra show 
clear power law behaviors: S{f) ~ f~f^. This behavior holds for all the other 
data. 

The non-stationary features have been investigated by varying the window-size 
on which the spectrum is calculated from 100 days up to the entire size of the 
time series. The power spectra coefficients (5 are calculated through a mean 
square regression in log-log scale. The values reported in Fig. 7 are the average 
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of the evaluated /3 over different windows and the error bars are their standard 
deviations. Fig. 7(a) refers to a time period between 1990 to 1996 whereas the 
Stock Market indices and Foreign Exchange rates (Fig. 7(b)) are analyzed over 
the time periods reported in Tables 1 and 2. Moreover, the averaged P values 
in a different time period, namely from 1997 to 2001 are reported in Table 7 
for Foreign Exchange rates. Stock Market indices and Treasury rates. These 
values differ from the spectral density exponent expected for a pure Brownian 
motion (/? = 2). However, we will shown in Section 8 that this method is 
biased and we indeed found power spectra exponents around 1.8 for random 
walks using three different random numbers generators. 

It must be noted that, the power spectrum is only a second order statistic 
and its slope is not enough to validate a particular scaling model: it gives only 
partial information about the statistics of the process. 

7.2 Scaling Spectral Density and Hurst Exponent 

For financial time series, as well as for many other stochastic processes, the 
spectral density S{f) is empirically found to scale with the frequency / as a 
power law: S{f) oc as already stated in the previous section. Here we use 
a simple argument to show how this scaling in the frequency domain should be 
related to the scaling in the time-domain. Indeed, it is known that the spec- 
trum S{f) of the signal X{t) can be conveniently calculated from the Fourier 
transform of the autocorrelation function (Wiener- Khinchin theorem). On the 
other hand, the autocorrelation function of X{t) is proportional to the second 
moment of the distribution of the increments which, from Eq. 2, is supposed 
to scale as K2 ~ r^^*-^-*. But, the components of the Fourier transform of a 
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function which behaves in the time-domain as are proportional to f~°'~^ in 
the frequency-domain. Therefore, we have that the power spectrum of a signal 
that scales as Eq. 2 must behave as: 

^(/) oc /-^^(^)-^ . (4) 



Consequently, the slope (3 of the power spectrum is related to the generalized 
Hurst exponent for g = 2 through: j3 = 1 + 2H{2). Note that Eq. 4 is obtained 
only assuming that the signal X{t) has a scaling behavior in accordance to 
Eq. 2 without making any hypothesis on the kind of underlying mechanism 
that might lead to such a scaling behavior. 

Wc here compare the behavior of the power spectra S{f) with the function 
j-2H{2)-i ^ff\^\(^ _ according to Eq. 4 - is the scaling behavior expected in the 
frequencies domain for a time series which scales in time with a generalized 
Hurst exponent H{2). We performed such a comparison for all the financial 
data and we report in Fig. 6 those for Stock Market indices for Thailand 
and JAPAN (in the time period 1997-2001). As one can see the agreement 
between the power spectra behavior and the prediction from the generalized 
Hurst analysis is very satisfactory. This result holds also for all the other 
data. Note that the values of 2H{2) + 1 do not in general coincide with the 
values of the power spectral exponents evaluated by means of the mean square 
regression. The method through the generalized Hurst exponent appears to be 
more powerful in catching the scaling behavior even in the frequency domain. 
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8 Monte Carlo Test of the Method 



In the hterature, the scahng analysis has been criticized for being biased. In 
order to test that our method is not biased we estimate the generahzed Hurst 
exponents for simulated random walks. We produce synthetic time series by 
using three different random number generators. We perform 100 simulations 
of random walks with the same number of data points as in our samples (991 
and 3118) and estimate the generalized Hurst exponents H{1) and H{2) and 
the power spectra exponents p. The results are reported in Table 5. 

In all the cases, H{1) and H{2) have values of 0.5 within the errors. Only when 
we consider uniformly distributed random numbers in the interval (0,1) (Rand 
which uses a lagged Fibonacci generator combined with a shift register random 
integer generator, based on the work of Marsagha [Marsagha and Zaman, 1994].) 
we obtain for H{1) of 0.47 ± 0.01, but also in this case H[2) is 0.5 within the 
errors. On one hand, this shows that our method is powerful and robust and is 
not biased as other methods are. On the other hand, the estimations of /3 from 
the power spectrum have values around 1.8 (instead of 2), showing therefore 
that this other method is affected by a certain bias. 

9 Conclusion 

By applying the same methodology to a wide variety of markets and instru- 
ments (89 in total), this study confirms that empirical scaling behaviors are 
rather universal across financial markets. By analyzing the scaling properties 
of the g'-order moments (Eq. 1) we show that the generalized Hurst exponent 
H{q) (Eq. 2) is a powerful tool to characterize and differentiate the struc- 
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ture of such scaling properties. Our study also confirms that qH{q) exhibits a 
non-linear dependence on q which is a clear signature of deviations from pure 
Brownian motion and other additive or uni-scaling models. 

The novelty of this work resides in the empirical analysis across a wide variety 
of stock indices that shows the sensitivity of the exponent H{2) to the degree 
of development of the market. At one end of the spectrum, we find: the Nasdaq 
100 (US), the ^&P500 (US), the Nikkei 225 (Japan), the Dow Jones Indus- 
trial Average (US), the CAC 40 (France) and the All Ordinaries index (Aus- 
tralia); all with H{2) < 0.5. Whereas, at the other end, we find the Russian 
AK&M, the Indonesian JSXC, the Peruvian LSEG, etc. (Fig. 4(b)); all with 
H{2) > 0.5. Moreover, we observe emerging structures in the scaling behaviors 
of interest rates and exchange rates that are related to specific conditions of 
the markets. For example, a strong deviation of the scaling exponent for the 3 
months maturity, which is strongly infiuenced by the central bank decisions. 
This sensitivity of the scaling exponents to the market conditions provides a 
new and simple way of empirically characterizing the development of finan- 
cial markets. Other methods usually used for controlling risk, like standard 
deviation or Sharpe Ratio are not able to provide such a good classification. 
The robustness of the present empirical approach is tested in several ways: by 
first comparing theoretical exponents with the results of Monte Carlo simula- 
tions using three distinct random generators, second by varying the maximum 
time-step (r^aa;) in the analysis, third by applying the Jackknife method to 
produce several samples, fourth by varying the time-window sizes to ana- 
lyze the temporal stability and fifth by computing results for detrended and 
non-detrending time series. We verify that the observed differentiation among 
different degrees of market development is clearly emerging well above the 
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numerical fluctuations. FinaUy, from the comparison between the empirical 
power spectra and the prediction from the scaling analysis (Eq. 4, Fig. 6) we 
show that the method through the generalized Hurst exponent describes well 
the scaling behavior even in the frequency domain. 
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Fig. 1. Kq{T) as a function of r in a log- log scale for the Nikkei 225 time series in 
the time period from 1990 to 2001 (r varies from 1 day to 19 days.). Each curve 
corresponds to different fixed values of q ranging from g = 1 to g = 3. In particular, 
the curve corresponding to q = \ (diamond markers), the one to q = 2 (square 
markers) and the curve to g = 3 (circle markers) are shown. 
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Fig. 2. The function qH{q) vs. q in the time period from 1997 to 2001. (a) JAPAN 
(Nikkei 225); (b) JAPAN (JPY/USD); (c) Thailand (Bangkok SET); (d) Thailand 
(THB/USD); (e) Treasury rates having maturity dates 9 = 10 years; (f) Eurodollar 
rates having maturity dates 6 = 1 year. For (f) the time period is 1990 - 1996. 
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Fig. 3. (a) The Hurst exponent H{1) for the Treasury and Eurodollar rates time 
series in the period from 1990 to 1996; (On the x-axis the corresponding maturities 
dates are reported.) (b) The Hurst exponent H(l) for the Stock Market indices and 
Foreign Exchange rates in the time period reported in Tabs. 1 and 2. (On the 
X-axis the corresponding data-sets are reported.) 
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Fig. 4. (a) The Hurst exponent H(2) for the Treasury and Eurodollar rates time 
series in the period from 1990 to 1996; (On the x-axis the corresponding maturities 
dates are reported.) (b) The Hurst exponent H(2) for the Stock Market indices and 
Foreign Exchange rates in the time period reported in Tabs. 1 and 2. (On the 
X-axis the corresponding data-sets are reported.) 
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Fig. 5. (a) The generalized Hurst exponent H{1) for the Stock Market indices in the 
whole time period (see Tab. 2) with its variation (black lines) obtained by using the 
Jackknife method and its variation (dashed lines) when time periods of 250 days are 
considered; (b) The generalized Hurst exponent H{2) for the Stock Market indices 
in the whole time period (see Tab. 2) with its variation (black lines) obtained by 
using the Jackknife method. The square points are the average values of H{1) and 
H{2) computed from a set of values corresponding to different Tmax- The error bars 
are their standard deviations. 
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Fig. 6. The power spectra of the Stock Market indices compared with the behaviour 
of f-^H{2)-i (straight hnes in log-log scale) computed using the Hurst exponents 
values in the time period 1997-2001; (a) Thailand (Bangkok SET) and (b) JAPAN 
(Nikkei 225). The line is the prediction from the generalized Hurst exponent H(2) 
(Eq. 4). 
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Fig. 7. (a) The averaged P values computed from the power spectra (mean square 
regression) of the Treasury and Eurodollar rates time series in the period from 1990 
to 1996; (On the x-axis the corresponding maturities dates are reported.) (b) The 
averaged /? values computed from the power spectra of the Stock Market indices and 
Foreign Exchange rates in the time period reported in Tabs. 1 and 2. The horizontal 
gray line corresponds to the value of /5 obtained from the simulated random walks 
reported in Table 5. (On the x-axis the corresponding data-sets are reported.) 
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Table 1 

Foreign Exchange rates (FX/USD). 
Country FX Time period 



Country FX Time period 



Hong Kong HKD 1990-2001 United Kingdom GBP 1990-2001 



Italy 



ITL 1993-2001 



France 



FRF 1993-2001 



Philippines PHP 1991-2001 Poland PLN 1993-2001 



Austraha AUD 1990-2001 



Peru 



PEN 1993-2001 



New Zealand NZD 1990-2001 Turkey TRL 1992-2001 



Israel ILS 1990-2001 



Thailand THE 1990-2001 



Canada CAD 1993-2001 



Mexico PESO 1993-2001 



Singapore SGD 1990-2001 Malaysia MYR 1990-2001 



Netherlands NLG 1993-2001 



India 



INR 1990-2001 



Japan JPY 1990-2001 Indonesia IDR 1991-2001 



Spain ESP 1990-2001 Taiwan TWD 1990-2001 



South Korea KRW 1990-2001 



Russia 



RUB 1993-2001 



Hungary HUF 1993-2001 Venezuela VEB 1993-2001 



Germany DEM 1990-2001 Brazil BRA 1993-2001 



Switzerland CHF 1993-2001 
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Table 2 

Stock Market indices (SM). 
Country 



SM 



Time period 



United States 



United States 



Japan 



Nasdaq 100 
S&P 500 
Nikkei 225 



Prance 
Australia 
United Kingdom 
Netherlands 

Germany 
Switzerland 
New Zealand 

Israel 
South Korea 
Canada 
Italy 
Spain 
Taiwan 



CAC 40 
All Ordinaries (AO) 
FTSE 100 
AEX 
DAX 
Swiss Market (SM) 
Top 30 Capital (T30C) 

Telaviv 25 (T25) 
Seoul Composite (SC) 
Toronto SE 100 (SE 100) 
BCI 30 
IBEX 35 
Taiwan Weighted (TW) 



1990-2001 



1987-2001 



1990-2001 



United States Dow Jones Industrial Average (DJIA) 1990-2001 



1993-2001 



1992-2001 



1990-2001 



1993-2001 



1990-2001 



1993-2001 



1992-2001 



1992-2001 



1990-2001 



1993-2001 



1993-2001 



1990-2001 



1990-2001 
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Table 2 (continued) 



Country 
Argentina 
Hong Kong 
India 
Brazil 
Mexico 
Singapore 
Hungary 
Poland 
Malaysia 
Thailand 
Philippines 



SM 



Merval (ME) 



Hang Seng (HS) 



Bovespa (BO) 



Mexico SE (MSE) 



Wig (WIG) 



Time period 



1993-2001 



1990-2001 



Bombay SE Sensex (BSES) 1990-2001 



1993-2001 



1993-2001 



All Singapore Shared (ASS) 1990-2001 



Budapest BUX (BUX) 1993-2001 



1991-2001 



KLSE Composite (KLSEC) 1990-2001 



Bangkok SET (BSET) 1990-2001 



Composite (CO) 



1990-2001 



Venezuela Indice de Cap. Bursatil (ICB) 1993-2001 



Peru Lima SE General (LSEG) 1993-2001 



Indonesia JSX Composite (JSXC) 1990-2001 



Russia AK&M Composite (AK&M) 1993-2001 
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Table 3 

Treasury rates {TRi{9)). 



1 3 months 7 7 years 

2 6 months 8 10 years 

3 1 year 9 30 years 

4 2 years 10 3 months (Bill) 

5 3 years 11 6 months (Bill) 

6 5 years 12 1 year (Bill) 
Table 4 

Eurodollar rates {ERi{0)). 

i e i e 

1 3 months 9 27 months 

2 6 months 10 30 months 

3 9 months 11 33 months 

4 12 months 12 36 months 

5 15 months 13 39 months 

6 18 months 14 42 months 

7 21 months 15 45 months 

8 24 months 16 48 months 
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Table 5 

Hurst exponents H{1) and H{2) and averaged /? values computed for random walks 
simulated by using three different random numbers generators: 1) Randn=Normally 
distributed random numbers with mean and variance 1; 2) Rand=Uniformly dis- 
tributed random numbers in the interval (0, 1) and 3) Normrnd=Random numbers 
from the normal distribution with mean and standard deviation 1. These are aver- 
age values on 100 simulations of random walks with 991 and 3118 numbers of data 
points. 

TV H{1) H{2) /3 



1) Randn 

991 0.50 ±0.01 0.50 ±0.01 1.8 ±0.1 

3118 0.50 ±0.01 0.50 ±0.01 1.80 ±0.03 

2) Rand 

991 0.47 ±0.01 0.49 ±0.01 1.8 ±0.1 

3118 0.47 ±0.01 0.50 ±0.01 1.80 ±0.03 
3) Normrnd 

991 0.49 ±0.01 0.49 ±0.01 1.8 ±0.1 

3118 0.50 ±0.01 0.50 ±0.01 1.80 ±0.03 
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Table 6 

Hurst exponents H{1) and H{2) for Foreign Exchange rates, Stock Market indices 
and Treasury rates in the time period from 1997 to 2001. 

Data H{1) H{2) Data H{1) H{2) 

Foreign Exchange rates 



HKD 



0.41 ±0.01 0.34 ±0.01 GBP 0.50 ± 0.02 0.48 ± 0.02 



ITL 



0.51 ±0.01 0.51 ±0.01 FRF 0.51 ± 0.01 0.51 ± 0.01 



PHP 



0.52 ±0.01 0.43 ±0.02 PLN 0.54 ± 0.01 0.50 ± 0.01 



AUD 



0.52 ±0.01 0.502 ±0.002 PEN 0.52 ± 0.01 0.41 ± 0.03 



NZD 



0.49 ±0.01 0.48 ±0.01 TRL 0.56 ± 0.01 0.44 ± 0.04 



ILS 
CAD 



0.48 ±0.02 0.47 ±0.02 THB 0.53 ± 0.01 0.50 ± 0.02 



0.51 ±0.01 0.48 ±0.01 PESO 0.53 ± 0.01 0.50 ± 0.01 



SGD 



0.50 ±0.01 0.47 ±0.03 MYR 0.51 ± 0.03 0.45 ± 0.05 



NLG 



0.51 ±0.01 0.51 ±0.01 INR 0.58 ± 0.02 0.53 ± 0.01 



JPY 
ESP 



0.50 ±0.01 0.49 ±0.01 IDR 0.56 ± 0.03 0.53 ± 0.03 



0.50 ±0.01 0.49 ±0.01 TWD 0.58 ± 0.01 0.51 ± 0.01 



KRW 



0.50 ±0.03 0.39 ±0.06 RUB 0.64 ±0.02 0.47 ±0.03 



HUF 



0.52 ±0.01 0.52 ±0.01 VEB 0.54 ±0.04 0.49 ± 0.02 



DEM 



0.51 ±0.01 0.51 ±0.01 BRA 0.59 ± 0.02 0.60 ± 0.01 



CHF 



0.51 ±0.01 0.50 ±0.01 
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Table 6 (continued) 



Data 



H{1) H{2) Data H{1) H{2) 



Stock Market indices 



Nasdaq 100 0.47 ± 0.01 0.45 ± 0.01 TW 0.53 ± 0.01 0.51 ± 0.01 



S&P 500 
Nikkei 225 
DJIA 
CAC 40 
AO 
FTSE 100 
AEX 
DAX 
SM 
T30C 
T25 
SC 
SE 100 
BCI 30 
IBEX 35 



0.47 ±0.02 0.44 ±0.01 ME 0.57 ± 0.01 0.53 ± 0.01 



0.46 ±0.01 0.43 ±0.01 HS 0.53 ± 0.01 0.49 ± 0.01 



0.49 ±0.01 0.464 ±0.004 BSES 0.54 ±0.01 0.52 ± 0.01 



0.47 ±0.02 0.46 ±0.02 BO 0.51 ± 0.01 0.48 ± 0.01 



0.49 ±0.02 0.46 ±0.03 MSE 0.57 ±0.01 0.52 ± 0.01 



0.46 ±0.02 0.44 ±0.01 ASS 0.57 ±0.01 0.54 ±0.02 



0.49 ±0.01 0.47 ±0.02 BUX 0.52 ± 0.01 0.49 ± 0.01 



0.50 ±0.01 0.47 ±0.01 WIG 0.49 ± 0.01 0.44 ±0.01 



0.50 ±0.02 0.48 ±0.02 KLSEC 0.60 ± 0.01 0.51 ± 0.02 



0.49 ±0.01 0.46 ±0.01 BSET 0.59 ± 0.01 0.55 ± 0.01 



0.53 ±0.01 0.51 ±0.01 CO 0.59 ± 0.01 0.54 ± 0.01 



0.53 ±0.01 0.51 ±0.01 ICB 0.61 ± 0.02 0.55 ± 0.02 



0.51 ±0.01 0.48 ±0.01 LSEG 0.61 ± 0.01 0.58 ± 0.01 



0.52 ±0.01 0.48 ±0.01 JSXC 0.57 ± 0.02 0.53 ± 0.02 



0.50 ±0.01 0.48 ±0.01 AK&M 0.65 ± 0.03 0.51 ± 0.01 
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Table 6 (continued) 



Data H{1) H{2) 
Treasury rates 

TRi 0.48 ± 0.01 0.44 ± 0.02 

TR2 0.55 ± 0.01 0.52 ± 0.02 

TRs 0.54 ± 0.01 0.52 ± 0.02 

Ti?4 0.53 ± 0.01 0.52 ± 0.02 

TR5 0.52 ± 0.01 0.50 ± 0.01 

Ti?6 0.51 ± 0.02 0.49 it 0.01 

TRy 0.49 ± 0.02 0.48 ± 0.01 

TRs 0.52 ± 0.01 0.50 ± 0.02 

TRg 0.51 ± 0.01 0.48 ± 0.01 

TRio 0.51 ± 0.01 0.48 ± 0.02 

Tiiii 0.56 ±0.01 0.54 ±0.02 

TR12 0.55 ± 0.01 0.53 ± 0.02 
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Table 7 

The averaged (3 values computed from the power spectra of the Foreign Exchange 
rates, Stock Market indices and Treasury rates in the time period from 1997 to 
2001. 

Data Averaged (3 Data Averaged (3 

Foreign Exchange rates 



HKD 



1.6 ± 0.2 GBP 1.79 ± 0.03 



ITL 



1.80 ±0.03 FRF 1.81 ±0.04 



PHP 



1.8 ±0.1 PLN 1.79 ±0.04 



AUD 



1.8 ±0.1 PEN 1.6 ±0.2 



NZD 



1.8 ±0.1 TRL 1.7 ±0.1 



ILS 
CAD 



1.8 ±0.1 THB 1.83 ±0.03 



1.80 ± 0.03 PESO 1.81 ± 0.04 



SGD 



1.81 ± 0.02 MYR 1.8 ± 0.1 



NLG 



1.81 ± 0.04 INR 1.8 ± 0.1 



JPY 



1.9 ± 0.1 IDR 1.83 ± 0.04 



ESP 



1.80 ± 0.04 TWD 1.8 ± 0.1 



KRW 



1.8 ±0.1 RUB 2.1 ±0.3 



HUF 



1.80 ± 0.03 VEB 1.8 ± 0.1 



DEM 



1.81 ±0.03 BRA 2.0 ±0.2 



CHF 



1.8 ±0.1 
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Table 7 (continued) 



Data 



Averaged P Data Averaged /3 



Stock Market indices 



Nasdaq 100 



1.7 ±0.1 TW 1.9 ±0.1 



S&P 500 
Nikkei 225 
DJIA 
CAC 40 
AO 
FTSE 100 
AEX 
DAX 
SM 
T30C 
T25 
SC 
SE 100 
BCI 30 
IBEX 35 



1.8 ±0.1 ME 1.8 ±0.1 



1.8 ±0.1 HS 1.8 ±0.1 



1.80 ± 0.03 BSES 1.82 ± 0.03 



1.8 ±0.1 BO 1.80 ±0.02 



1.8 ±0.1 MSE 1.9 ±0.1 



1.81 ±0.03 ASS 1.9 ±0.1 



1.8 ±0.1 BUX 1.82 ±0.04 



1.8 ±0.1 WIG 1.8 ±0.1 



1.8 ±0.1 KLSEC 1.8 ±0.1 



1.8 ±0.1 BSET 1.9 ±0.1 



1.9 ±0.1 CO 2.0 ±0.2 



1.9 ±0.1 ICB 2.0 ±0.2 



1.9 ±0.1 LSEG 2.0 ±0.2 



1.9 ±0.1 JSXC 1.9 ±0.1 



1.8 ± 0.1 AK&M 1.9 ± 0.2 
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Table 7 (continued) 

Data Averaged /3 Data Averaged (3 
Treasury rates 

TRi 1.8 ±0.1 TRy 1.9 ±0.1 

TR2 1.83 ±0.04 TRs 1.9 ±0.1 

TR3 1.86 ±0.05 TRg 1.8 ±0.1 

ri?4 1.88 ±0.06 TRio 1.82 ±0.04 

TR5 1.9 ±0.1 TRu 1.85 ±0.04 

TRq 1.9 ±0.1 TR12 1.9 ±0.1 
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